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Abstract. Article describes a class of efficient algorithms for 3SAT. 



Introduction 

This paper presents the compatibility matrix method for 3SAT. The compatibil- 
ity matrix is an encoding of 3SAT instances using the contradictions between their 
clauses. This encoding is the first major result of this work. 

The compatibility matrix encoding allows parallel testing of all guesses. This 
paper describes several algorithms which perform the parallel testing of all guesses 
in polynomial time indeed. The algorithms transform the compatibility matrix into 
a similar encoding of the satisfying true-assignments. These efficient algorithms for 
3SAT are the second major result of this work. 

Let us describe 3SAT as follows. Suppose, / is a given conjunctive normal form 
(CNF) with clauses of length three or less. Let 

(0.1) / = c 1 Ac 2 A...Ac m 

In formula /, clauses Ci are disjunctions of three or less non-complimentary literals 
(different Boolean variables or their negations): 

c» € {aV/3V7, aVft a | a,/3,7 e | j = 1,2, . . . , n}}, 

- where bj are Boolean variables. Suppose, we have to decide whether there does 
exist such a true-assignment which would satisfy formula /: 

true = f(bi, b 2 , ■ ■ ■ , b n )\ bj=T . e {f a i se j rue y 

Such a decision problem is a 3SAT instance. 

Cook proved that 3SAT is a NP-complete problem pQ. So, the whole history 
of the NP-completeness theory is the history of 3SAT. As a matter of fact, 3SAT 
was one of the first four NP-complete problems discovered pQ . One might say that 
3SAT is the most beautiful NP-complete problem discovered. 



1. Compatibility matrix 

Let Ti be the truth-table for clause a, i = 1, 2, . . . , m. Let's arbitrarily enumerate 
the true-assignments for the arguments of Ci and let's write Tj in the following form: 



Ti 


True- assignments 


Ci 


1 

2 


1- st true-assignment 

2- nd true-assignment 


Value of Ci on the 1-st true-assignment 
Value of Ci on the 2-nd true-assignment 
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Let Si^ be u-th string from truth-table T. By definition, any two strings Si^ 
and Sjv are compatible if the following two conditions are satisfied: 

Condition 1: The /i-th value of q is true and the z^-th value of Cj is true; 

Condition 2: The /i-th true-assignment in table T and the ^-th true- assign- 
ment in table Tj do not contradict each other. In other words, if clauses 
Ci and Cj share a variable, then this variable has to have the same true- 
assignment in both strings and Sj v . 

When at least one of these conditions is not satisfied, then the strings are incom- 
patible. Let us emphasize, i and j may be equal in this definition, i.e. strings Si U 
and Sj V may be from the same truth-table. 

Let's build a compatibility box for each clause couple (cj, Cj), i,j = l,2,...,m. 
Let ki and kj be the lengths of clauses Ci and Cj, appropriately. Then, the compat- 
ibility box for clauses ci and Cj is a Boolean matrisQ dj — {x liu ) 2ki y 2 k j with the 
following elements: 



true. Strings and Sj V are compatible 
false, Otherwise 



Having all m compatibility boxes built, we can aggregate them in a box matrix C: 

C ( C% j ) rn X m 

- where the size is given in "boxes" . Matrix C is a Boolean box matrix. For 3SAT, 
size of C is 

rn m 

(^2 fel ) x (J2 2kz ) < ( 8 ™) x (8m) = 0(m x m), 

1=1 i=l 

- where ki is length of clause Ci. We call matrix C a compatibility matrix for formula 

itm 

Let us notice that the compatibility matrix is a symmetric matrix: 



(1.1) C = C T , dj = C. 



Also, let us notice that, due to Condition 2 for the string compatibility, diagonal 
boxes Cn are diagonal matrices, i.e. all non-diagonal elements in these boxes arc 
false: 



(i.2) a. 



( t\ false ■ ■ ■ \ 
false t2 

v '■■ ••• / 



, tfj_ € {true, false} 



Because clauses in formula 10.11 are disjunctions, there is the only one value false 
among values t^ in matrix [L^l 

As usual, we can use a (0, Inversion of the compatibility matrix. In the version 
of matrix C, values true are replaced with 1, and values false are replaced with 0. 



^Boolean matrix is a matrix whose elements are true or false. 
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2. Solution grid and depletion 

Traditionally, guessing in 3SAT is the checking whether a true-assignment sat- 
isfies formula lO.il But, any true-assignment for the variables in formula 10.11 con- 
tains/consists of the true- assignments for clauses Cj. So, any true-assignment pro- 
duces in the compatibility matrix a grid of elements, one element per compatibility 
box: 

7 = {aVij I M = v = i,j = 1,2, . . . ,m}, 
- where x^ij is the (/i, ^)-element of compatibility box Cy , and 7 is a grid of such 
elements in compatibility matrix C, one element per compatibility box. 

Based on the grid presentation, we can state that any true-assignment satisfying 
formula 10.11 is presented in compatibility matrix C with a grid of elements, one 
element per compatibility box, which all are true. And visa versa, if there is such 
a grid of elements in the compatibility matrix, one element per compatibility box, 
whose all elements are true, then formula lO.ll is satisfiable. We call such a grid 
which consists of all irite-elementH, one element per compatibility box, a solution 
grid. 

Based on property 11.21 of the compatibility matrix, we could even remove the 
"one-element-per-box" requirement from the solution grid definition. 

Theorem 2.1. Formula \0.1\ is satisfiable iff compatibility matrix for that formula 
contains a solution grid. 

Solution grid is a pattern of satisfiability in the compatibility matrix encoding 
of formula 10.11 And theorem 12.11 reduces the guessing to the searching the matrix 
for the pattern. In the search, any heuristic is welcome. 

Our method is a detection of those elements of the compatibility matrix which 
do not belong to any solution grid and the assigning to them value false or 
depending on the presentation of the compatibility matrix. We call the replacement 
of a true-element with value false a depletion of the compatibility matrix. 

Theorem 2.2. Suppose, there is a filter which finds in the compatibility matrix 
for formula \0.1\ all such "true "-elements which do not belong to any solution grid. 
Suppose, all these elements were replaced with value "false". Then, formula \0.1\ is 
satisfiable iff there is not any compatibility box filled with "false" entirely. 

In fact, for unsatisfiable formula 10. \\ the matrix emerging from the depletion 
described in theorem 12.21 will be filled with false entirely. 

A compatibility box entirely filled with false we call a pattern of unsatisfiability. 
One might say that our method is a search of the compatibility matrix for the 
pattern. For the method, it does not matter what particular filter is used for the 
depletion. For the filter, the only objective is to preserve at least one solution grid 
when the grids exist. Then, no information will be lost in the sense of computational 
complexity. 

3. Basic algorithm 
The following algorithm is a depletion filter. We call it a basic algorithm: 



True-element is an element which is equal true, and /a/se-element is an element which is 
equal false. 
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Init: Build the compatibility matrix for formula lO.il 

Let s = 1. Let C s be the compatibility matrix. Let C S ij be the compatibility 

boxes, where i,j = 1,2, ... , m. 
Depletion: Let s = s+ 1. Calculate compatibility matrix C s as follow^: 

m 

(3.1) C S ij = A Cs-i^kCs-i^kj , 

fe=l 

- where C s jj is the (i, j)-th compatibility box of C s , i, j = 1, 2, . . . , m. 
Let 

Xsfivij be the (/i^)-th element of compatibility box C S ij. Then, due 
to formula [3. 11 the clement is 



(3-2) ^sfiuij — ^\ ( \J % s — l,[iaik A X s — \_QLukj ) ■ 

fc=l l<a<2 3 

Due to propertv ll.il of the compatibility matrix, the last expression can be 
rewritten in any of the following ways: 

(3.3) ^s^iuij j\(\J •Es — l,ai±ki A X s — l.avkj*) j\(\J -^s—l.^iaik ^ ^ s — l^uajk) • 

k a k a 

Iterations: If C s ^ C s _i, then go to the previous step. Otherwise, continue. 

Decision: If C s is entirely filled with false, then formula [0.11 is unsatisfiable 
- decision "NO". Otherwise, formula [0.11 is satisfiable - decision "YES". 
Obviously, when during an iteration the pattern of unsatisfiability arises (value 
false fills a compatibility box entirely), then the algorithm can be stopped with 
decision "NO" . Because, formula 13.11 will propagate this value false all over the 
compatibility matrix during the next two iterations, at most. 

For 3SAT, the computational complexity of the algorithm can be estimated as 
0(m 5 ): there is 0(m 2 ) iterations (each iteration depletes at least one element 
from the compatibility matrix); there is m 2 depletions on each iteration; and each 
depletion takes time 0(m). 

Theorem 3.1. Any "true "-element in the matrix emerging after running of the 
basic algorithm belongs to a solution grid. Formula \0.1\ is unsatisfiable iff the basic 
algorithm produces the pattern of unsatisfiability for this formula. 

Proof. Clauses in formula 10.11 can be expressed in their disjunctive forms (DF). 
Depending on their length, the clauses are 

a = false V a 

, .s a V (3 = false V a f\ (3 V a A [3 V a /\ (3 

^ ' ' aV/3V7= false V aA/3A7 V aA/3A7 V aA(3Aj\/ 

V aA/3A7 V a A (3 A 7 V a A (3 A 7 V a A (3 A 7 

- where a, /3, and 7 are literals. 

Let's call all conjunctions in DF 13.41 the triads (value false in the formulas is 



Here, product of two Boolean matrices A = (a^v) and B = (b^) of the appropriate sizes (the 
number of columns in A has to be equal to the number of rows in B) is the following Boolean 
matrix: 

AB= (\/a ik Ab kj ). 
k 

Conjunction of Boolean matrices of the same size is the matrix of conjunctions of the appropriate 
elements of the matrices. 
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a triad, too). For each clause, there is a natural one-to-one relation between the 
clause's triads and true-assignments: that true-assignment which makes a clause 
equal false relates to triad false; the rest of the true-assignments are in the relation 
with those triads to which they deliver value true. Let's enumerate the triads in 
DF 13.41 using that one-to-one relations: triad is that triad which relates to /i-th 
true-assignment for clause Cj. 

Let's replace clauses in formula rO.ll bv their DF I3.4| and let's open all parentheses 
in the formula. The result of that is the following DF: 

m 

/= V A r ^ 

/il,/i2,---,/im t=l 

- where the disjunction is taken over all appropriate combinations of indexes fii. 
Then, formula [0.1 1 is satisfiable iff there is at least one satisfiable conjunction 

m 

(3.5) T(fj,x,fi2,---,H m ) = /\T fli i 

i=l 

In other words, formula 10.11 is satisfiable iff there is at least one conjunction 13.51 
without complimentary literals. 

Now, let's write the following symbolic box matrix Y: 

Y — (^ij)mxm, 

- where Y^ are the boxes. The symbolic boxes in matrix Y are as follows: 

Yij — (r^j A T U j) 2k . x2 kj , 

- where fej and kj are lengths of clauses Ci and Cj appropriately. 

Let's calculate powers of matrix Y using formulas 13.11 and 13.21 

rn 

Y = Y, Y = ( /\ Y Si ikY s ,kj)mxm, s > 1 
k=l 

- where boxes Y^,*,* are the appropriate boxes in box matrix Y s . Due to formula 
13.21 elements of matrix Y s = (y a i3)mxm are some DF of the following kind: 

(3.6) y a p = \J 7) Vi A t^ 2 ,a 2 A... 

CTl ,CT2 .... 

- where (ui, (72, ■ ■ •) are all appropriate combinations of the triads' indexes. 

The most important for us is to notice that, starting from some power s = 0(m 2 ) 
of matrix Y, each conjunction in DF 13.61 contains at least m different triads. But, 
any conjunction of more than m different triads contains at least two different triads 
from the same clause. But, any such conjunction equals false because there are 
complimentary literals in those different triads from the same clause. Thus, starting 
from that power, all DF 13.61 get finalized as follows: 

(3.7) y al3 = \J Tim, H2, . ■ ■ , Hm) 

- where the disjunction is taken over all appropriate combinations of indexes /Lt,. 

Let's rewrite r(/xi, fj,2, ■ ■ ■ , Mm) as follows: 

m 

r(jUl,/X 2 , . . . ,Mm) = f\ T Mii = 
i=l 
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= (v,iAt Wj i) A (t^.iAt^^) A ... A (^.lAr^m) A 

A (t^jAvj) A (t M2 , 2 At M2 , 2 ) A ... A (t M2i 2 A T^ m , m ) A 

A ... A 

A (r Mm , m A r Mli i) A (r^ m , m A r M2i2 ) A ... A(r Mim A^ B , m ) 

Each parenthesis in the last presentation is an element of the compatibility matrix 
for formula fO.il 

A T u j — X fj,t/ij 

- where x^uij is the (p,v)-the element in compatibility box dj. Element x^ V ij is 
true iff triads T^i and r U j contain no complimentary literals. Thus, conjunctions [375] 
are associated with the grids of the compatibility matrix's elements, one element 
per compatibility box: 

2^2, Ml, 2,1 A 2^2^2,2,2 A ... A X M2:Atmi 2,m A 

A 

%[i m ,fii,m,l A Xfj !rn _fj !2 ^ rn _2 A ... A X^ m : fj, m ^m.m 

Among all conjunctions 13.51 only those are equal true whose elements constitute a 
solution grid in the compatibility matrix: 

, {J>2 : ■ • ■ : /^m) — *£jUi ,/A2sl)l — •EfAi ,/^2 , 1,2 — ■■• — ^jUl ,l,m 



r(/Ui,^2, • ■ ■ ,Mm) 



(3.8) 



A 
A 
A 



true 



□ 



Compatibility matrix is an encoding of formula lO.il Theorem 13 . 1 1 shows that the 
basic algorithm transforms this encoding into an encoding of those true-assignments 
which satisfy the formula. Each true-element in the resulting matrix belongs to a 
solution grid. The true-assignments appropriate to any resulting true-element are 
a part of a satisfying true- assignment. These parts can be used to compute the 
whole satisfying true-assignments by the method of self-reducibility. 

Let us notice that, due to formula ^. 21 the length of conjunctions in DF 13. 61 grows 
exponentially over power of matrix Y. So, the expected number of iterations in the 
basic algorithm can be estimated as O(logm). Thus, the expected running time of 
the algorithm can be estimated as 0(m 3 log to). 

Also, in the case of unsatisfiable formula 10.11 number to in all estimations of 
time can be replaced by the number of clauses in the minimal unsatisfiable sub- 
formula of formula 10.11 That is due to the conjunctions in the formula 13.11 the 
basic algorithm restricted to those compatibility boxes which are built for a mini- 
mal unsatisfiable sub-formula will deplete less true-elements but it will produce the 
pattern of unsatisfiability. 

Before going further, let us notice that, due to Condition 1 of the strings com- 
patibility, the true-assignments, which deliver to their clauses value false, can be 
missed when building the compatibility matrix, because they will produce in their 
compatibility boxes the rows and columns entirely filled with false. So, the com- 
patibility matrix can be built solely on the contradictions between true-assignments 
satisfying different clauses separately. 
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4. Asynchronous depletion 

Obviously, depletions 13.11 could use the current values of the compatibility boxes 
which were computed during the current iteration. Such an asynchronous depletion 
can be expressed with the following back-feeds: 

(4.1) C i3 = C i3 A C ik Ckj, 1 < i, j, k < m 

Really, diagonal boxes Cu in the compatibility matrix are the diagonal Boolean 
matrices, anqj 

Cij = CaCij = CijCjj, Ca < CijCij T , Cjj < C^Cij 

Thus, depletions 13.11 are conjunctions of depletions 14. II 
Then, the basic algorithm can be reformulated as follows: 

Init: Compose the compatibility boxes for formula lO.ll and select a schema 
to iterate index-triplets 

(4.2) {i,k,j): l<i,j,k<m 

The schemes can be context-sensitive and have other specifics. For example, 
a schema could "visit" the triplets with different frequencies or even miss 
some of them. The simplest schemes just cycle trough all index-triplets; 

Depletions: Follow the selected iteration schema and deplete the appropriate 
compatibility boxes using transformations 14.11 until each of these boxes 
became an invariant under the depletions. 

Decision: If the pattern of unsatisfiability emerges during the depletions, 
then formula [0.11 is unsatisfiable. Otherwise, the formula is satisfiable. 

The freedom of selection the iteration schema creates a family of depletion filters: 
the asynchronous algorithms differ one from another by their iteration schemes. 
Theorem 13.11 implies 

Theorem 4.1. Let an iteration schema traverse all index-triplets \JJ^ Then, for- 
mula \0.1\ is unsatisfiable iff depletions \4-l\ produce the pattern of unsatisfiability for 
this formula. 

For 3SAT, computational complexity of the asynchronous algorithms can be es- 
timated as 0(m 3 ): there are ra 3 index-triplets l4~2l during iterations 14. 1[ each of 
these triplets involves O(l) Boolean operations to deplete the appropriate compat- 
ibility boxes of size 2 3 x 2 3 or less. 

Obviously, depletions 14.11 will detect the pattern of unsatisfiability (if there is 
one) in time 0(s 3 ), where s is the number of clauses in the shortest unsatisfiable 
sub- formulas of formula lOTTl 

5. Triangular algorithms 

Due to the symmetry of the compatibility matrix, any asynchronous depletion 
can be restricted to those compatibility boxes which are located over/under the 
major diagonal of the matrix. For example, the compatibility boxes allocated over 
the major diagonal of the compatibility matrix are described with the following 
index-triplets: 

(5.1) 1 < i < k < j <m 



false = false, false < true, true = true. 
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For these index-triplets, depletions 14. II have to be rewritten as follows: 

!Cij — C'ij A C'ikCkj 
Cm = Cife A CijC^j 
Ckj = Cfcj A C ik C i: j 

Theorem 14.11 implies the following. 

Theorem 5.1. Let an iteration schema traverse all index-triplets fXil Then, for- 
mula Un\ is unsatisfiable iff devletions \5.£\ vroduce the pattern of unsatisfiability for 
this formula. 

Any such algorithm we call a triangular algorithm. 

6. Square algorithms 

For the (0, Inversion of the compatibility matrix, depletions 13. 31 can be rewritten 
"sub-tropically" : 



(6.1) x^ij = min{x^ij, min{l, 



- ^maxfi^i^j | 1 < a < 2 fe "} 



m 

(3=1 



}}, 



- where x^ij is the (/x, z^)-th element of compatibility box Cy, hp is the number of 
literals in clause c@, and [*J is the whole part of a number. 

Obviously, all the above algorithms do work with depletions 16.11 and have the 
same computational complexity. The above theorems imply 



Theorem 6.1. Formula \ 0. 1\ is unsatisfiable iff a square algorithm produces the 
pattern of unsatisfiability for this formula. 

We call such algorithms the square algorithms because of the following reason. 

Obviously, the compatibility matrix for formula 10.11 contains a solution grid iff 
there are m such permutation matrices X%, X2, ■ ■ ., and X m which would satisfy 
the following system of inequalitie^]: 



(6.2) XidjXf > 



( 1 


V = 



i,j = 1,2,..., m 



- where ki and kj are the numbers of literals in clauses Ci and Cj appropriately. 
This system just reflect that fact that, for a satisfiable formula [0.11 there is such a 
re-enumeration of true-assignments in truth-tables Tj which would move a solution 
grid's elements in the upper left corners of their compatibility boxes. 

System [6.21 is a system of integral bilinear inequalities. Iterations 16. II iust solve 
this system. Formula 16.11 iteratively compares the squares of the compatibility 
matrix with the squares of that matrix on the right side of system RT21 and depletes 
the compatibility matrix appropriately. 

The result of such depletions is a general solution of system 16.21 in the form of 
the compatibility matrix encoding of all true-assignments satisfying formula 10.11 



5 



For matrices A = (a^) and B = (b^) of the same size, relation A > B means that 

> 1 
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Conclusion 

This paper described on high-level several polynomial time algorithms for 3SAT. 
The algorithms are based on a original encoding of 3SAT instances - a compati- 
bility matrix. The encoding reflects the contradictions between true-assignments 
satisfying different clauses separately. Such a "hologram" of 3SAT instances allows 
parallel testing of all guesses. The algorithms transform the compatibility matrix 
into a matrix encoding of the satisfying true-assignments. 

Basically, the compatibility matrix presents 3SAT instances in the form a jigsaw 
puzzle, and the algorithms simulate the human activities during the solving of the 
puzzle. 

For some of the algorithms, there are demos at [2]. Readers can use these demos 
to probe the algorithms with their own 3SAT instances. 
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